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, On the basis of the classical theory of envelope, we formulate the renormalization 

group (RG) method for global analysis, recently proposed by Goldenfeld et al. It is 
clarified why the RG equation improves things. 
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1 Introduction 



a: 

Recently, Goldenfeld et al|lj have proposed a new method based on the renormalization group 
(RG) equation 0, |J to get the asymptotic behavior of solutions of differential equations. 
The method is simple and has a wide variety of applications including singular and reductive 
perturbation problems in a unified way. However, the reason is obscure why the RG equation 
can be relevant and useful for global analysis: The RG equation is usually related with the 
scale invariance of the system under consideration. The equations which can be treated in 
the RG method are not confined to those with scale invariance 0. Actually, what the RG 
method does in [|lj] may be said to construct an approximate but global solution from the ones 
with a local nature which were obtained in the perturbation theory; the RG equation is used 
to improve the global behavior of the local solutions. This fact suggests that the RG method 
can be formulated in a purely mathematical way without recourse to the concept of the RG. 
A purpose of the paper is to show that this is the case, thereby reveal the mathematical 
structure of the method. 

Our formulation is based on the classical theory of envelopes @. As everybody knows, 
the envelope of a family of curves or surfaces has usually an improved global nature compared 
with the curves or surfaces in the family. So it is natural that the theory of envelopes may 
have some power for global analysis. One will recognize that the powerfulness of the RG 
equation in global analysis and also in the quantum field theory [pi [3l is due to the fact that 
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it is essentially an envelope equation. We shall also give a proof as to why the RG equation 
can give a globally improved solution to differential equations. 

In the next section, a short review is given on the classical theory of envelopes, the notion 
of which is essential for the understanding of the present paper. In §3, we formulate the RG 
method in the context of the theory of envelopes and give a foundation to the method. In 
§4, we show a couple of other examples to apply our formulation. The last section is devoted 
to a brief summary and concluding remarks. 



2 A short review of the classical theory of envelopes 

To make the discussions in the following sections clear, we here give a brief review of the 
theory of envelopes. Although the theory can be formulated in higher dimensions^, we take 
here only the one-dimensional envelopes, i.e., curves, for simplicity. 

Let {C T } T be a family of curves parametrized by r in the x-y plane; here C T is represented 
by the equation 

F(x,y,T) = Q. (2.1) 
We suppose that {C T } T has the envelope E, which is represented by the equation 

G(x,y) = 0. (2.2) 
The problem is to obtain G(x, y) from F(x, y, r). 

Now let E and a curve C TQ have the common tangent line at (x, y) = (xq, yo), i.e., (xq, yo) 
is the point of tangency. Then xq and yo are functions of r ; xo = 4>(to), yo = V ; ( r o), and of 
course G(xo,yo) = 0. Conversely, for each point (xo,yo) on E, there exists a parameter r . 
So we can reduce the problem to get r as a function of (x , yo); then G(x, y) is obtained as 
F(x,y,T(x,y)) = G(x,y)B t (x ,?/o) can be obtained as follows. 

The tangent line of E at (x ,yo) is given by 

f(r )(x - x ) - (j>'{ro){y - yo) = 0, (2.3) 

while the tangent line of C T0 at the same point reads 

F x (x , y , t )(x - x ) + F y (x , yo, r )(y - y ) = 0. (2.4) 

Here F x = dF/dx and F y = dF/dy. Since both equations must give the same line, 

F x {x , y , to)4>'{to) + F y (x , y , To)ip'(To) = 0. (2.5) 

1 Since there is a relation G(xo,yo) = between xo and yo, To is actually a function of xq or yo. 
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On the other hand, differentiating F(x(tq), y(r ), r ) = with respect to r , one has 

F x (x ,yo,r )(j)'(T ) + F y (x , y , r )^'(r ) + F To (x , y , r ) = 0, (2.6) 

hence 

j-i / x _ dF(x ,y ,r ) 

F T0 [xo, |/o, r ) = 5 = 0. (2.7) 

One can thus eliminate the parameter r to get a relation between x and y , 

G(x,y)=F(x,y,r (x,y)) = 0, (2.8) 
with the replacement (x ,y ) — > (x,y). G(x,y) is called the discriminant of F(x,y,t). 

Comments are in order here: (i) When the family of curves is given by the function 
y = f(x,r), the condition Eq.(2.7) is reduced to df/dr = 0; the envelope is given by 
y = f(x, t (x)). (ii) The equation G(x, y) — may give not only the envelope E but also a set 
of singularities of the curves {C T } T . This is because the condition that dF/dx = dF/dy = 
is also compatible with Eq. (2.7). 

As an example, let 

y = f{x, r) = e— (1 - e • (x - r)) + e~ x . (2.9) 
Note that y is unbound for x — r — > oo due to the secular term. 

The envelope E of the curves C T is obtained as follows: From df /dr = 0, one has r — x. 
That is, the parameter in this case is the x-coordinate of the point of the tangency of E and 
C T . Thus the envelope is found to be 

y = f(x,x)=e- e >-e-*. (2.10) 

One can see that the envelope is bound even for x — > oo. In short, we have obtained a 
function as the envelope with a better global nature from functions which are bound only 
locally. 

As an illustration, we show in Fig.l some of the curves given by y — f(x,r Q ) together 
with the envelope. 



Fig.l 



3 Formulation of the RG method based on the theory 
of envelopes 

In this section, we formulate and give a foundation of the RG method [0] in the context of the 
classical theory of envelopes sketched in the previous section. Our formulation also includes 
an improvement of the prescription. 



3 



Although the RG method can be applied to both (non-linear) ordinary and partial dif- 
ferential equations, let us take the following simplest example to show our formulation: 

d doc n , > 

He + e n + z = °- <w 

where e is supposed to be small. The solution to Eq.(3.1) reads 



x{t) = Aexp(-~t) sin(y l-jt + 6), (3.2) 
where A and 9 are constant to be determined by an initial condition. 

Now, let us blindly try to get the solution in the perturbation theory, expanding x as 

x(t)=x (t) + exi(t) + e 2 x 2 (t) + (3.3) 
where x n (n = 0, 1, 2...) satisfy 

x + x = 0, x n+ i + x n+1 = -x n . (3.4) 

Thus xo = A Q sm(t + 9o), and then x\ + x± — —A cos(t + 6*o), and so on. Then we get for 
xi and special solutions 

xi(t) = (t-t )sm(t + e ), 

x 2 (t) = ^{(t-t ) 2 sm(t + 9 )-(t-t )cos(t + 9 )}. (3.5) 

o 

Here we have intentionally omitted the unperturbed solution from x n (t) (n = 1,2, ...). Al- 
though this prescription is not adopted in the proceeding calculations are simplified with 
this prescription; see also §4.i It should be noted that the secular terms have appeared in the 
higher order terms, which are absent in the exact solution and invalidates the perturbation 
theory for t far away from £ . 

Inserting Eq.(3.5) into Eq.(3.3), we have 

x(Mo) = A sm(t + 9 ) -ey(t-t o )sin(t + o ) 

+e 2 ^{(t - t ) 2 sin(t + 9 ) - (t - t ) cos(t + 9 )} + 0(e 3 ) (3.6) 



Now we have a family of curves {Ct }t given by functions {x(t,to)}t parametrized with 
to- They are all solutions of Eq. (3.1) up to 0(e 3 ), but only valid locally, i.e., for t near to- 
Let us find a function x E {t) representing the envelope E of {C to }t - 

2 It is amusing to see that the unperturbed solution in the higher order terms x n (n = 1,2, ...) is anal- 
ogous to the "dangerous" term in the Bogoliubov's sense in the quantum-filed theory of superfluidity and 
superconductivity^) . 
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According to the previous section, we only have to eliminate to from 

0, (3.7) 



dx(t,t ) 



dt 

and insert the resultant t (t) into x(t,t ). Then we identify as x E (t) = x(t,t (t)). It will be 
shown that x E (t) satisfies the original differential equation Eq. (3.1) uniformly for Vt up to 
0(e 4 ); see below. 

Eq.(3.7) is in the same form as the RG equation, hence the name of the RG method[|l|]. 
In our formulation, this is a condition for constructing the envelope. 

Here comes another crucial point of the method. We assume that A and B are func- 
tionally dependent on t ; 

A = A (t ), e = e (t ), (3.8) 

accordingly x(t,t ) = x(t, A (t ), 9o(t ), t ). Then it will be found that Eq. (3.7) gives a 
complicated equation involving A (t ), 9 (t ) and their derivatives as well as t . It turns out, 
however, that one can actually greatly reduce the complexity of the equation by assuming 
that the parameter t coincides with the point of tangency, that is , 

t = t , (3.9) 

because A (to) and 6*o(to) can be determined so that to = t. We remark here that the 
meaning of setting to = t is not clearly explained in [Ij], while in our case, the setting has 
the clear meaning to choose the point of tangency at t = to-§ 

From Eq.'s (3.7) and (3.9), we have 

d-Ao . d6 e 2 

-d^ + eAo = °> ^ + 8 = °- (3 ' 10) 
Solving the simple equations, we have 

M*o) = Ae~^l\ e (t Q ) = - E \ + 6, (3.11) 

o 

where A and 9 are constant numbers. Thus we get 

e e 2 

x E (t) = x(t, t) = Aexp(--t) sin((l - -)t + 9). (3.12) 



Noting that yl — e 2 /4= 1 — e 2 /8 + 0(e 4 ), one finds that the resultant envelope function 
x E (t) is an approximate but global solution to Eq.(3.1); see Eq. (3.2). In short, the solution 
obtained in the perturbation theory with the local nature has been "improved" by the 
envelope equation to become a global solution. 



3 It is interesting that the procedure to get the envelope of x(t, Ao(to), Oo(to), to) assuming a functional 
dependence of Aq and 6*o on to is similar to the standard prescription in which the general solution of a 
partial differential equation of first order is constructed from the complete solution. H] 
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There is another version of the RG method]!]], which involves a "renormalization" of the 
parameters. We shift the parameter for the local curves as follows: Let r be close to t, and 
write t — to = t — r + r — t . Then putting that 



A(t) = A Q (t Q )Z(t ,T), Z(t ,r)=l--(T-t ) + -(r-to 

Z o 

0(r) = 6 (t ) + 56, 59 = -^(r-t ), (3.13) 

o 

we have 

x(t,r) = A(r)sin(t + 9(T))-e^p-(t-T)sm(t + 9(T)) 

+e 2 ^l {(t - rf sin(t + 6{t)) - (t - r) cos(t + 9{r))} + 0(e 3 ), (3.14) 

o 

where 

x(t,t) = A(r)sm(T + 9(r)). (3.15) 

Then the envelope of the curves given by {x(t, r)} r will be found to be the same as given in 
Eq. (3.12). 

This may concludes the account of our formulation of the RG method based on the 
classical theory of envelopes. However, there is a problem left: Does x E (t) = x(t,t) indeed 
satisfy the original differential equation? In our simple example, the result Eq.(3.12) shows 
that it does. It is also the case for all the resultant solutions worked out here and in 
[H]. We are, however, not aware of a general proof available to show that the envelope 
function should satisfy the differential equation (uniformly) up to the same order as the 
local solutions do locally. We give here a proof for that for a wide class of linear and non- 
linear ordinal differential equations (ODE). The proof can be easily generalized to partial 
differential equations (PDE). 

Let us assume that the differential equation under consideration can be converted to the 
following coupled equation of first order. 

^ = F(q(t),t;e), (3.16) 

where 'q = (qi, q%i ■ ■ •) and F are column vectors. It should be noted that F may be a 
non-linear function of q and t, although in our example, 

« = * = *■ F = (_ 9l _ f92 ). < 3 - 17 > 

i.e., F is linear in q. We also assume that we have an approximate local solution q(t, to) 
around t = to up to 0(e n ); 

^ = F(q,t;e) +0(e"). (3.18) 
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One can see for our example to satisfy this using Eq.(3.6). 
The envelope equation implies 

<9q(t, t ) 







at to = t. With this condition, q E (t) corresponding to x E (t) is defined by 

q B (t)=q(t, t). 



(3.19) 



(3.20) 



It is now easy to show that q B (t) satisfies Eq.(3.16) up to the same order as q(t, to) does: 
In fact, for Vt 



dq E {t) 



dt 



dq(t, t ) 



t=t 



dt 



dq(t,t ) 



t=t 



dt 



dq(t,t ) 



t=t 



dt 



(3.21; 



t=t 



where Eq.(3.19) has been used. And noting that F(q B (to), t ; e) = F(q(t , t ), t ; e), we see 
for Vt 



dt 



F(q E (t),t;e) + 0(e n ) 



(3.22) 



on account of Eq. (3.18). This completes the proof. It should be stressed that Eq.(3.22) is 
valid uniformly for Vt in contrast to Eq.(3.18) which is valid only locally around t = t . 



4 Examples 



Let us take a couple of examples to apply our formulation. These can be converted to 
equations in the form given in Eq. (3.16). 



4.1 A boundary- layer problem 



The first example is a typical boundary-layer problemQ: 

d2 y i n i \ dy i n 
dx z dx 



(4.1) 



with the boundary condition y(0) = 0, y(l) = 1. The exact solution to this problem is 
readily found to be 



y(x) 



exp(— x) — exp(— x/e) 
exp(— 1) — exp(— 1/e) 



(4.2) 
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Now let us solve the problem in the perturbation theory. Introducing the inner variable 
X by eX — x ||, and putting Y(X) = y(x), the equation is converted to the following; 

dxi + dx= ~<Jx + Y) - (4 ' 3) 

Expanding Y in the power series of e as Y = Y + eY\ + e 2 Y 2 + one has 

Yo+Yo = 0, 

y;' + y; = -Fo'-Fo, 

:■ (4-4) 

Here, Y' = dY/dX etc. To solve the equation, we set a boundary condition to Y(X) and 
Y (X) at X = X ; 

Y(X) = Y (X ) = A , (4.5) 
where Xq is an arbitrary constant and A$ is supposed to be a function of Xq. 

For this problem, we shall follow the prescription given in fl[ for the higher order terms. 
Then the solutions to these equations may be written as 

Y (X) = A - B e- {x - Xo \ 

Y 1 (X) = -A (X-X )-(B + C )(e-( x - x ^-l). (4.6) 
Defining A = A + e(B + C ) and B = B + e(B + C ), we have 

Y(X, X Q ) = A - Be~ ix - Xo) - eA{X - X ) + 0(e 2 ). (4.7) 
In terms of the original coordinate, 

y(x, x Q ) = Y(X, X ) = A - Be- {x - Xo)/e - A(x - x ) + 0(e 2 ), (4.8) 
with xo = X /e. 

Now let us obtain the envelope Y E (X) of the family of functions {Y(X, X )}x each of 
which has the common tangent with Y E (X) at X = X . According to the standard procedure 
to obtain the envelope, we first solve the equation, 

dY 

0, with X = X, (4.9) 



dX 

and then identify as Y(X,X) = Y E (X). 
Eq. (4.9) claims 

A' + eA = 0, B' + eB = 0, (4.10) 

with the solutions A(X) = Aexp(-eX), B(X) = B exp(-X), where A and B are constant. 
Thus one finds 

Y E (X) = Y(X,X) = A(X) - B(X) = Ae~ tX - Be~ x . (4.11) 
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In terms of the original variable x, 

y B (x) = Y B (X) = Aexp(-x) - 5exp(~). (4.12) 

It is remarkable that the resultant y E {x) can admit both the inner and outer boundary 
conditions simultaneously; y(0) = 1, y(l) = 1. In fact, with the boundary conditions we 
have A = B = l/(exp(— 1) — exp(— 1/e)), hence y E (x) coincides with the exact solution y(x) 
given in Eq. (4.2). 

In Fig. 2, we show the exact solution y(x) and the local solutions y(x,xo) for several x$: 
One can clearly see that the exact solution is the envelope of the curves given by {y(x, x )} xo . 



Fig.2 



A comment is in order here: If we adopted the prescription given in §3 for the higher 
order terms, the perturbed solution Yi(X) reads Yi(X) = —A (X — X ); note the boundary 
condition Eq.(4.5). Then the proceeding calculations after Eq. (4.6) would be slightly 
simplified. 



4.2 A non-linear oscillator 

In this subsection, we consider the following Rayleigh equation]?], [jj, 

y + y = e(y-^f). (4.13) 



Applying the perturbation theory with the expansion y = y + eyi + e 2 y2 + ■ ■ •, one has 

R R 3 

y(t,t ) = R sm(t + 6o) + e{(^ - -£-)(t - t ) sm(t + 6 ) 

+^(cos3(t + ^o)} + 0(6 2 ). (4.14) 

Here we have not included the terms proportional to the unperturbed solution in the higher 
order terms in accordance with the prescription given in §3, so that the following calculation 
is somewhat simplified than in Jl]. Furthermore, the result with this prescription will coincide 
with the one given in the Krylov-Bogoliubov-Mitropolsky method]?]], as we will see in Eq. 
(4.18). 

Eq. (4.14) gives a family of curves {Ct }t parametrized with to- The envelope E of 
{C to }t with the point of tangency at t = to can be obtained as follows: 

^|M = , (4 . 15 ) 
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with to = t. Assuming that Rq and 9q are ~ 0(e) at most, we have 



2 

the solution of which reads 



A, = - 9 = 0, (4.16) 



run = _f° (4.17) 

'exp(-et) + Rl{l - exp(-et))/4 



with Rq = -Ro(O) and 9q = constant. Thus the envelope is given by 

y E (t) = y(t,t) = R (t)sm(t + e ) + e^^(cos3(t + 6 )) + O(e 2 ). (4.18) 

This is an approximate but global solution to Eq. (4.13) with a limit cycle in accordance 
with the result given in [^]. We note that since Eq.(4.13) can be rewritten in the form of 
Eq.(3.16), Eq.(4.18) satisfies Eq.(4.13) up to 0(e 2 ). 



5 A brief summary and concluding remarks 



We have given a geometrical formulation of the RG method for global analysis recently 
proposed by Goldenfeld et al[l]: We have shown that the RG equation can be interpreted 
as an envelope equation, and given a purely mathematical foundation to the method. We 
have also given a proof that the envelope function satisfies the same differential equation up 
to the same order as the functions representing the local curves do. 

It is important that a geometrical meaning of the RG equation even in a generic sense 
has been clarified in the present work. The RG equation appears in various fields in physics. 
For example, let us take a model in the quantum field theory ||; 

C = \{d,ct>f - ±<f + at., (5.1) 

where c.t. stands for counter terms. The true vacuum in the quantum field theory is 
determined by the minimum of so called the effective potential V(0 C )|| In the one-loop 
approximation, the renormalized effective potential reads 

V(0 C , M) = + ^4(ln -4 - -), (5.2) 

where M 2 is the renormalization point. To see a correspondence to the envelope theory, one 
may parametrize as <p 2 = exp t and M 2 = exp i > then one sees that hi(f) 2 /M 2 becomes a 
secular term t — t . In the quantum field theory, one applies the RG equation to improve 



the effective potential as follows 0, |10] 

dV 



()if2 0, with M 2 = 4>1 (5.3) 
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One sees that this is the envelope equation! The resultant "improved" effective potential is 
found to be 



V impr (0 c ) = V(0 C , <j> c ) = if^jr- (5-4) 



16tt 2 



>rl) 



Thus one can now understand that the "improved" effective potential is nothing but the 
envelope of the effective potential in the perturbation theory. One also sees the reason why 
the RG equation with <p c = M can "improve" the effective potential. Then what is the 
physical significance of the envelope function Vi mpr (0 c )? One can readily show that for VM 



dV((f>, 



c j impr 



<9V(0 C , M) 



(5.5) 

<t> c =M 



owing to the envelope condition Eq. (5.3). This implies, for example, that the vacuum 
condensate (fi c that is given by dVi mpr /d(f>1 = is correct up to the same order of ^-expansion 
in which the original effective potential is calculated; this is irrespective of how large is the 
resultant C . Detailed discussions of the application of the envelope theory to the quantum 
field theory will be reported elsewhere [|TT|] B 



The RG equation has also a remarkable success in statistical physics especially in the 
critical phenomena ||. One may also note that there is another successful theory of the 
critical phenomena called coherent anomaly method (CAM) [^3| . The relation between CAM 
and the RG equation theory is not known. Interestingly enough, CAM utilizes envelopes 
of susceptibilities and other thermodynamical quantities as a function of temperature. It 
might be possible to give a definite relation between CAM and the RG theory because the 
RG equation can be interpreted as an envelope equation, as shown in this work. 

Mathematically, it is most important to give a rigorous proof for the RG method in 
general situations and to clarify what types of differential equations can be analyzed in this 
method, although we have given a simple proof for a class of ODE's. We note that the proof 
can be generalized to partial differential equations, especially of first order with respect to a 
variable [O]. One should be also able to estimate the accuracy of the envelope theory for a 
given equation. We hope that this paper may stimulate studies for a deeper understanding 
of global analysis based on the theory of envelopes. 
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Figure Captions 

Fig.l A family of functions and its envelope: 

The thin lines show y = exp(— er )(l — e(x — r )) +exp(— x) with r = 0.2, 0.4, 0.4, 0.6 
and 0.8, which are attached to the respective lines. The thick line shows the envelope 
y = exp(— ex) — exp(— x). (e = 0.8.) 

Fig. 2 The thin lines show y(x, xo) = A(xo) — B{xq) exp(— (x — x )/e) — A{xq){x — x ) with 
xq = 0.2, 0.4 and 0.8, which are attached to the respective lines. The thick line shows 
y(x) given in Eq.(4.2). (e = 0.1.) 
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